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Abstract 



The recent discovery of the B c meson by the CDF collaboration and proposed 
new experiments at Fermilab and CERN motivate new theoretical studies of 
the B c system. Here we investigate the radiative leptonic decay B c — > jlu. 
This process is an important background to the annihilation process B c — > Iv, 
which will be used to extract the B c decay constant. We perform a model- 
independent calculation, based on QCD, of the partial width and various 
kinematic distributions. We also examine the decay within the framework of 
NRQCD, an effective field theory of nonrelativistic quarks, generalizing the 
NRQCD Lagrangian to include external sources for the weak and electromag- 
netic currents. Finally, we will show how NRQCD reproduces the correct 
position of the B* pole in the emission of very soft photons. 
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I. INTRODUCTION 

The B c is probably the final long lived pseudoscalar meson which will be found in our 
lifetimes. Its recent discovery at CDF M, and the future prospect of thorough experimental 
studies at LHC and BTeV, motivates extensive theoretical studies of this system 0. In 
particular, the B c offers a unique possibility to study the effects of weak interactions in a 
quarkonium-like environment. Since it is composed of heavy quarks of two different flavors, 
the B c is stable against strong annihilation decays and its dynamics approaches a simple 
perturbative limit as m c , mj — > oo. In this limit, the B c is a very compact bound state of 
a c and b quark, with a small admixture of non-perturbative higher Fock states containing 
gauge bosons and light quark-antiquark pairs. This admixture is small because soft nonper- 
turbative gluons, with large Compton wavelengths, have little overlap with the compact B c 
state. It is useful to study the B c system, like quarkonium, in the framework of an effective 
nonrelativistic quantum field theory 0. 

The advantages of an effective field theory approach are both conceptual and quantita- 
tive. While one would intuit that the m c , rrib — ► oo limit is one in which a description based 
on a constituent quark model should work reasonably well, an effective field theory based 
on the operator product expansion puts this intuition on a rigorous basis. It is unsurprising, 
of course, that the results we obtain in this paper are similar in structure to those of the 
quark model, since the quark model does respect the symmetries of the heavy quark limit. 
What an effective field theory does, that the quark model does not, is provide an oganized 
expansion within which the leading corrections to this limit may be accomodated. While 
we will not compute higher order corrections in this paper, we will use the power counting 
of the effective field theory to estimate their size, thereby casting light on the accuracy of 
our results. 

One basic characteristic of the B c meson is its decay constant, defined by 

(n\h»i 5 c\B c (jp)) = if B y. (i.i) 

This is a QCD definition, although in quark models and in the nonrelativistic limit f Bc may 
be identified with the value of the wave function at the origin. The decay constant probes 
the strong QCD dynamics which is responsible for the binding of the quark-antiquark state. 
The most straightforward way to determine f Bc would be to measure the purely leptonic 
partial widths, such as B c — ► lv\, where I = e, fi,r ||. The rate for this process is given by 

ns c - iu) = ^|y c6 | 2 /i<4- f i - 4-) 

8tt c m Bc V m B J 

However, the practical usefulness of this method is limited, because the decaying meson 
is spinless and this mode is helicity suppressed. A helicity flip on an external lepton line 
is required, leading to a suppression of the rate by an additional factor m 2 /m 2 Bc . This 
suppression is 6 x 10~ 8 for B c — > ev e , and 3 x for B c — > fiu^. The only leptonic mode 
with a substantial branching fraction is B c — > tu t , but this is difficult to observe because 
the r must be reconstructed from its decay products. 

The helicity suppression can be overcome if there is third particle in the final state. In 
particular, adding a photon does not change the fact that the decay rate is proportional 

1 



(1.2) 



to the decay constant. A naive estimate suggests that for I = fi the additional electro- 
magnetic coupling is effectively compensated by the lifting of the helicity suppression, since 
(a/An)/(m Bc /m 2 ^) ~ 2. For I = e, the same counting suggests that the radiative leptonic 
mode dominates the purely leptonic decay. If branching ratios as small as that for B c — > fiu^ 
are eventually measured, then there are two consequences. First, B c — > "ffiu^ will be an im- 
portant background (or tool) for the extraction of fs c - Second, B c — > ^ev e will be observable 
even though B c — > ev e is not. For both of these reasons, it is important to understand the 
radiative decay process. 

Radiative leptonic B c decays already have been studied using quark potential models 
and QCD sum rules 0. We will comment on these approaches later and relate them to 
our own, which will be based on QCD and its nonrelativistic expansion, NRQCD. The 
corresponding decay for systems with one heavy and one light quark, B — > also has 
been examined 0, employing insofar as possible a heavy fermion expansion for the b quark. 
Unfortunately, that decay is dominated by photons radiated from the u, the effect of which 
is impossible to compute model-independently. The advantage of the B c system is that it 
can be treated systematically in an expansion in v, the nonrelativistic three-velocity of the b 
and the c in their mutual bound state. The expansion may be used to justify the application 
of perturbative QCD to this decay [[J. In addition, we note that, in contrast to B — > lu, •ylu, 
the decays B c — > lu, •ylu are not CKM-suppressed. 

The value of NRQCD is that it provides a rigorous counting of powers of v, which 
may be applied to expectation values of heavy quark operators in external states QQ 
where ttiq, mqi ^> Aqcd- This power counting is related to the heavy quark effective theory 
(HQET) derivative expansion in powers of D^/tuq, although it differs in some details. What 
NRQCD and HQET share is a "spin symmetry" , by which the magnetic interactions of the 
heavy quarks decouple as l/niq. In general, the power of v which NRQCD assigns to 
an operator depends on the matrix element which is being taken. Therefore the NRQCD 
expansion must be performed on the matrix elements, not just on the operators. By contrast, 
the HQET power counting is identical for all bound states to which it is applied, which are 
of the form Qq, where m q < Aqcd- Although the HQET expansion is technically also an 
expansion of matrix elements, there is therefore no confusion in thinking of it equally as 
applying to the operators themselves. 

In fact, we will be able to treat our application of NRQCD to the B c system as an operator 
expansion as well, because of a number of simplifications which apply to the process under 
consideration. First, we will work only to leading order in the NRQCD expansion, with no 
dependence on the \bcg) higher Fock components for which the power counting in v can be 
subtle (and, of course, interesting). Second, the only external states which will be important 
are the B c and the B*, both of which are dominated by the same spatial configuration of 
the be pair. Therfore, the relevant leading matrix elaments are related by heavy quark 
spin symmetry. The only matrix element which will appear is the one which defines fs c in 
NRQCD, namely 

im B J Bo = (0\xli>c\B c ) , (1.3) 

where Xb and ip c are the two-component NRQCD field operators. Third, the final state 
consists of no particles which are strongly interacting. As a result, the operator product 
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FIG. 1. Feynman diagrams for be — > 'ylvi 

expansion, which NRQCD helps to organize in the case of hadronic decays, is here a trivial 
affair. 

We will first study the decay B c — > Iwy in perturbative QCD, which is fairly simple 
since there are no infrared divergences or other subtleties. We will then match the QCD 
result onto NRQCD and identify the leading operator which contributes to the annihilation 
of the B c . The expansion turns out to be straightforward, but it is useful to see how the 
parton model result emerges in the nonrelativistic limit. Finally, we will show how NRQCD 
reproduces the correct position of the B* pole in the emission of very soft photons. 

II. THE DECAY IN PERTURBATIVE QCD 

We begin by calculating the rate for B c — > Iwy perturbatively in full QCD. Since there 
are no strongly interacting particles in the final state, the QCD operator product expansion 
is trivial. The matrix element for B c — > Iwy is governed by the matrix element for bc — > lw~f, 
followed by the projection of the be onto the B c state. At leading order, we take the quarks 
to be in the (leading) 1 5'o Fock configuration. 

The computation is simple, involving the set of diagrams presented in Fig. [T| The decay 
amplitude can be written as 



AQ>c[ x S Q ] -> Iwy) 



Gp - 
^ C 6-^e*w(p l/ )7 M (l - 7 5 )v(pi) (0| &7 a (l - 7s)c \bc) 



(2.1) 



p c k p b k pik 



e Qc T nva _ eQbPc 



Pck 



Pbk 



e QlPc pew/j 

Pik 



where 



(2.2) 
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Here e M and k^ are the polarization and momentum of the photon. We now anticipate the 
results of the next section, in which we will verify that the current 670,(1 — 75)0 is the leading 
contribution to the decay, to justify the evaluation of the B c matrix element flLip . Defining 
the leptonic matrix element = uj^l — 75) u, we find 



A(B C -> luj) = ~iV cb ^ fB c m B J,e*J^v^ + 



eQcPc _ eQbPb 
p c k p b k 



+ B^kvg^ - kN v ) - B 2 ie^v a k^ , (2.3) 

where t> M is the four-velocity of the B c , B\ = eQ c /2p c k — eQb/2pbk + eQi/2pik, and B 2 = 
eQ c /2p c k + eQb/2pt,k + eQi/2pik. This expression simplifies somewhat if we choose the 
second gauge condition e • v — 0, which we can do since the photon is on shell. 

This is the leading result in the nonrelativistic expansion. It differs from existing cal- 
culations in two respects |5|||. First, it is model-independent. Second, radiation from the 
lepton leg is included, which is required for gauge invariance of the amplitude. At this point, 
it is straightforward to derive the total rate and the photon and lepton energy distributions. 
Clearly, they depend only on the one parameter fg c . Since the same parameter enters the 
expression for the rate of B c — > lis, we will normalize the rate to that for B c — > lis, and 
keep the dependence on the fermion charges explicit so that the various contributions can 
be examined. Summing over the photon and lepton polarizations and integrating over the 
phase space, we find the decay rate 



T{B C -> lisj) 
T(B C 



lis) 



a 


m B c ~ 




Air 


mf _ 


m k ( 



9ml 



IQl 2 Q b Qi 



9 m\ 



9 rribm r 9 m c m b, 



2 QcQi _j_ 2 Qf 



9m 2 Bc 



(2.4) 



We see that there is no interference between the photon emitted from the charm and bottom 
legs. This is a consequence of the anticorrelation of the spins of the two quarks in the 
pseudoscalar B c . The dimensionless coefficient r\ = (a/7i)(m 2 Bc /mf) is r M = 2.1 for / = fi 
and r e = 8.8 x 10 5 for / = e. Taking nib = 4.8 GeV, m c = 1.5 GeV and ms c = 6.3 GeV, we 
find 



T{B C -> lisj) 
T(B C -> lis) 



(0.19 Ql + 1.96 Ql - 0.29 Q b Qi + 0.93 Q C Q X + 0.22 Q 



n (0.02 + 0.87 - 0.10 - 0.62 + 0.22) 

0.40 n. 



(2.5) 



Hence the rate for B c — > fiisj is approximately 80% of that for B c — > fiis, while B c —>■ ewj 
dominates B c — > eis. The doubly differential spectrum in x = 2E 1 /rriB c and y = 2Ei/rriB c is 



1 dr 

r dxdy 



18(1 



- x)x~ 2 {/ C {1- 6/i 6 ) + 2^(2 - Q/2 C + 9/4)) 1 - 6/j b ) (1 — y) 2 

+ 4/^(1 - 3/i c )(l - x - y) 2 + 9^ 2 C (2 - 2x + x 2 - Ay + 2xy + 2y 2 ] , (2.6) 



where fib = mb/m-Bc and /i c = m c /mB c - The normalized photon energy spectrum follows 
the simple shape 
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FIG. 2. Photon energy spectrum in B c — > Ivj decay. 

i dr . , 

-- = 6.(1 - z) . (2.7) 

The photon energy spectrum is shown in Fig. |||. Note that there is no soft divergence in 
the limit E 7 — > 0. This is because the helicity suppression for B c decay can only be lifted 
by recoil effects or a heavy quark spin flip from a magnetic transition, both of which vanish 
with E T The lepton energy spectrum, whose analytic form is 



p j- = 9(/i c 2 (l - 6fi b ) + 2fi 2 b (2 - 6/i c + 9^ 

2/i c 2 y(l - 6/i 6 )(l -y)+ n 2 b y(24 + 72^ - 63/$, - 20y - 72^ c + 60/i c y) (2.8) 
- 2(1 - y) ln(l - y) (2^((6/x c - 2)(3 - y) - 9^(2 - y)) + ^( 6 /i 6 - 1)(1 - y) N 1 



is shown in Fig. [5| Unlike the photon energy spectrum, the lepton energy spectrum depends 
on the mass ratios /t c and //&. Finally, we compute the branching ratio for B c — > Ivy. 
Taking for the B c lifetime the CDF central value t(B c ) = 0.46 ps PJ, and using the estimate 
f Bc = 375 MeV [|, we obtain 

B(B C -»■ ei/7) ~ B(B C -> /iz/7) ~ 4.4 x 10" 5 . (2.9) 

It is certainly a challenge to observe such rare processes, even with the large B c samples 
which one hopes will be available at future experiments such as BTeV and LHC 

This result is similar to what one would find using a constituent quark model. However, in 
a number of respects it goes beyond the quark model framework. First, the result holds even 
when the photon is hard and the virtual c or b quark is far from its mass shell. Second, the 
it is rigorously true in the limit m c , — > 00. Third, it is possible to improve systematically 
the accuracy of this result by including both QCD radiative corrections and higher dimension 
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FIG. 3. Lepton energy spectrum in B c — > Zz^7 decay. 

operators in NRQCD. The QCD calculation provides insight into both the success and the 
limitations of the quark model picture. 

Note that for very soft photons, the shape of the spectrum is dominated by the nearest 
be resonance, the B*, and the parton-hadron duality which underlies the operator product 
expansion breaks down. (Since the splitting between the B c and the B* is caused by a 
hyperfine interaction, it is much smaller than the splittings between the B c and all other 
B** excitations.) This corner of phase space will be examined in more detail in Section IV. 
At larger photon energies, all virtual B** states contribute equally to the shape, but the 
integration over the momenta of the leptons smears the effect of the resonances into a smooth 
result. By the usual application of global parton-hadron duality, the smeared spectrum 
should be reproduced by perturbative QCD (supplemented by a hierarchy of nonperturbative 
contributions from operators of higher dimension). Note that the kinematics forbids the 
production of on-shell B** intermediate states. Nor do on-shell cc bound states contribute 
to this exclusive process via the diagram in Fig. 0a, since the physical photon has k 2 = 0. 

III. THE NRQCD EXPANSION FOR THE B c 

We now turn to the nonrelativistic limit of the QCD answer, matching onto a tower of 
NRQCD operators. The purpose is to identify the leading contribution to the B c matrix 
elements which contribute to B c — > ^Iv. (The result, which turns out to be simple, was 
already anticipated in the previous section.) We will write our expansion in terms of effec- 
tive fields ^(x), which are related to the usual QCD fields Q(x) by a Foldy-Wouthuysen 
transformation, 

Q(x) = exp{Hp ± /2m Q } x fy(x) = exp{i~D ■ ^ /2m Q } x #(x) , (3.1) 
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where = — v^v ■ D. In the rest frame of the B c and in the Dirac representation, \I/ 
may be decomposed as 



(3.2) 



where ip and x are t ne two-component quark and antiquark fields, which the Foldy- 
Wouthuysen transformation disentangles to order 1/rriQ. The field ip annihilates quark 
states, while x creates antiquark states. Note that the projection operators P± = (1 ±^)/2, 
where v M = (1, 0, 0, 0) is the four-velocity of the B c , project out separately the ip and x parts 
of It is often convenient to use this covariant (and representation- independent) form of 
the projection. 



The transformation (|3.1| ) can be used to rewrite QCD operators in terms of the effective 
fields This is one step in the matching of QCD onto NRQCD. For example consider the 
expansion of the weak current bT^c = 67^(1 — 7 5 )c to order l/m c fc. In the four-component 
notation, we find 



(3.3) 



In the two-component notation, this becomes 



67^(1 



7 5 )c 



fi() 



1 



- — xl lD a Wc + - — Xb° %L ] Wc 



x\iD ip c -\ 

2rrib 2m c 



X h iD 3 ip c 



(3.4) 



—ie jkl x\iD k a l i) c 



which is somewhat more cumbersome but also more explicit. These leading terms in the 
expansion in \jm^ c are also the leading terms in the velocity expansion; the first omitted 
terms are of order v 2 . A more complete calculation would include radiative corrections as 
well, in which case the various operators which appear on the right hand sides of Eqs. ( p.3|) 
and ( p.4| ) would develop a dependence on the renormalization scale. Note that the covariant 
derivative includes both gluon and photon gauge fields. 

The effective fields \l/ realize explicitly the heavy quark spin symmetry which emerges as 
the magnetic interactions of Q decouple in the limit itlq — > 00. An immediate implication of 
this symmetry in NRQCD is that the B c and the B* are degenerate. Another is the equality 
of their decay constants, 



im B J Bc = (0\xl^Pc\B c ) = ^(OlxhVcl^) + 0(v 2 



(3.5) 



A particularly nice realization of this symmetry is available in terms of the four- component 
effective fields. One can assemble the B c and the B* into a single "superfield" 



H 



(3.6) 
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Matrix elements which conserve the spin symmetry then take the form 



(0|^r* e |BW>= % -f Bc m B l^[TH Bo ], 

(sW|^ c r* c |BW) = -m B Ti [H B TH Bc 
(B^\%T^ b \B^) = m B Tr [h b TH Bc ] . (3.7) 

The leading contribution to the mass difference between the B c and the B* comes from the 
chromomagnetic operator 

2 = -^-^ c a^G a ^ c + -^-%a a(3 G a ^ b . (3.8) 
4m c 4m b 

The mass difference A = m B * — m Bc may be written in terms of the matrix elements of 2 , 

A = (B*\0 2 \B*) - (B c \0 2 \B c ). (3.9) 
This relation will be useful later, when we consider the effect of the B* pole. 

A. NRQCD calculation and matching 

Due to the simplicity of the final state, we can match the currents of QCD directly onto 
NRQCD, rather than the decay rates. For example, Lorentz invariance restricts the form of 
the NRQCD current which annihilates the B c to 

(0\f\B c ) = -g^ IdiOlx&cW + -^-(0\xl D 2 ^ C \B C ) + ...), (3.10) 



m 



red 



where m rec i = m b m c / {m b + m c ) is the reduced mass of the be pair. At leading order, this 
reduces to the matrix element of a single universal NRQCD operator which can be either fixed 
by other experimental measurements or related to the B c wave function at the origin jnj . 

We would like to generalize the effective Lagrangian of NRQCD to include the effects of 
the weak interactions. This can be achieved by introducing external sources. The effective 
NRQCD Lagrangian can be separated into two sectors representing flavor conserving (jCfcw) 
and flavor changing {Cpch) interactions, as well as a source part £ 5 , 

£-NRQCD = £-FCon J r C-FCII + C-S ■ (3-H) 

The sources represent external perturbatively interacting fields. In addition to the electro- 
magnetic field, in this calculation it is convenient to treat the leptonic current as an external 
source field responsible for the flavor changing interactions. Each part of this Lagrangian 
can be organized by the velocity expansion. At leading order, the flavor conserving part can 
be written 

£ F Con = Xl {i A - D 2 /2m a ) Xa + £ (i A + D 2 /2m a ) ^ a 

■ v 2 ip a - Xa° ■ $2Xo) , (3- 12) 



On 



2m n 
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where the covariant derivatives include external scalar s, D t = dt + a$s and vector vx, 
D = d + ia^vi sources, and summation over quark flavors a = {c, b} is understood. Here 
af represent the coefficients of NRQCD operators, which are normalized such that af = 
1 + 0(a s ). Note that in this notation, the field x creates an antiquark, so if the Lagrangian 
( |3.12j ) were normal ordered then the kinetic term for the antiquark field would have the 
opposite sign. At leading order, the flavor changing part of Eq. (|3.11|) reads 



*-FCh 



(3.13) 



where Si and V\ represent flavor-changing external sources. It is important to remember 
that the power counting rules for the Lagrangian with external sources are different from the 
usual NRQCD Lagrangian, in that one counts only the powers of quark fields. Therefore, 
the next-to-leading Lagrangian contains covariant derivative insertions into the leading order 
Lagrangian, 



r>NL 

L 'FCon 



bf_ 

8m;? 



5mf 



Xct 



D x 



a 



Xc 



x a) ip, 



vi. (3.14) 



For the external electromagnetic field one identifies s\ and V\ with the external elec- 
tromagnetic potentials, while v 2 and v% represent external magnetic and electric fields, 
v\ = eQaW = \eQ a e ijk F jk and v\ = eQo0 = eQ a F 0i . At next-to-leading order, the 
flavor changing Lagrangian is 



L 'FCh 



iD 



\ 



, cr tp p S 2 - —Xo 
771/3 I 2 



iD iD 1 
+ 

m a nip 



I'M 



+ b ie iM xl[ — ~ — \*%Vl+ h.c. 
2 \ m n rap 



(3.15) 



We will identify Si and Vi with the time and space components of the leptonic current 
= V tA (G F /V2)u'y li (l - 7 5 )w = V cb (G F /V2)l^. As defined in Eqs. (CT) and (p!5|) , the 

coefficients satisfy bf = 1 + 0(a s ). 

A set of NRQCD Feynman rules follows directly from this Lagrangian and can be used 

to calculate the set of diagrams of Fig. [p. The calculation of diagram Fig. |TJc is identical to 

the one in QCD and yields 



Ac = a 3 V c 



cb 



G F eQi 
2V2( Pl k) 



{2v»rf + {kv)sT ~ Wv v - ie^v a kp] emO\xl^ c \bcCSo)), (3.16) 



The amplitudes of diagrams Fig. [TJa and Fig. |I]b can be computed using Feynman rules 
derived from Eqs. (|3.12|) and Q3.13|) , 



Gi 



A a , b = -ta 2 a,V c b^e ijk rk^* k (0\xl^ c \bcCS )) 



eQb 



+ 



eQ c 



2m b {pi - Pi/2m b ) 2m c [j>% - p 2 2 /2m c ) 



(3.17) 
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where p\ = E b + E 17 p x = p b + k, pi = E c - E 1 , p 2 = p c - k. 

To compare Eq. (|3.17|) to the corresponding expression in full QCD, it is convenient to 
adopt a covariant notation and choose the gauged e • v = 0. Then it can be rewritten as 



A ■ T / G p 1 

A a ,b = ^a 2 a A V cb -j= — 



eQ b eQ c 



2m b 2m r 



^VXM0|x^cM%)>. (3.18) 



In addition, there is also a local interaction coming from a 64 term of the flavor changing 
part of the Lagrangian ( |3.15|) , 



^i + £i\^ c \bcCS )) (3.19) 



m b m c 



This contribution involves operators with covariant derivatives. Recalling that covariant 
derivatives contain electromagnetic source fields whose couplings can be treated perturba- 
tively, we obtain 

A s = -h,V cb % - \sT - vW) etl,(0\xl^\bcCSo))- (3.20) 

2V2 \m c m b J 

Note that this term has no dependence on the photon energy. Adding the contributions 
( p. 161 ), ( 3.18|) and (|3.20|) and comparing with the QCD result ( |2.1|) , we reproduce at leading 



order the anticipated matching conditions bf — af — 1. 

We hasten to add a remark here. At leading order, the only contribution that survives in 
the denominators of Eq. ( |3.17| ) is E^ = kv, matching the propagator of full QCD. However, 
in NRQCD this calculation is only good for _E 7 m c ^ b , where the heavy quarks are close 
to their mass shell. There is a substantial portion of the final state phase space in which 
the propagating heavy quark is relativistic and therefore not described by the formalism of 
NRQCD at the leading order. In a formal sense, the sum of an infinite number of NRQCD 
operators is required to fully describe the propagation of the relativistic quark. We can 
separate the region in which NRQCD is valid from the one in which it is not by introducing 
a factorization scale /1 -C m Cjb . Then the amplitude ( |3.18| ) is generated by the NRQCD 
Feynman diagrams in Fig. |l| only for E 1 < \i. For the region E 1 > //, we must add a local 
"operator" to the Lagrangian, 

- + ^) Xli; (3.21) 



a/2 kv \2m b 2m b 

to reproduce NRQCD amplitude in this kinematic region. This object is an operator in the 
NRQCD quark fields, but a c-number with respect to the external electromagnetic sources. 
Note that its inclusion is required by the matching conditions of QCD onto NRQCD, to 
account for the situation in which the intermediate heavy quark is far from its mass shell . 



Clearly, the matching procedure can be performed in any gauge; in NRQCD, this gauge choice 
is most convenient and quite common. 
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This is entirely in accord with what one usually finds upon applying the operator product 
expansion, that highly virtual intermediate states generate operators of higher dimension in 
the low energy theory. Note that since the this new object is to be included only in the case 
of large photon energy E 1 > fi, it could equally well be written in terms of a local operator 
built out of the photon field strength F^. We have chosen to write Eq. ( |3.21| ) in this form 



so that it may be compared transparently to the amplitudes which contribute at small E 7 . 

It is now clear that the leading contribution to the decay B c — > '-flu indeed comes from 
the dimension three matrix element (O|x]/0c|-E?c) = irriB c fB c , as anticipated in the previous 
section. This entirely unsurprising result has now been placed on firm theoretical foot- 
ing, as the leading term in a systematic expansion. Furthermore, we are in a position to 
estimate the size of the most important corrections. From the expansion (|3.10| ), we see 
that there is a relativistic correction to (O|xj/0 C |.B C ) of relative order v 2 by NRQCD power 
counting ||. This is a contribution to /# c , however, which cancels in the ratio ( |2.5|) . The 
leading correction which is particular to B c — > ^Iv comes from the dimension five operator 
(Olxls"? ' B^d-Bc) /ml, generated by attaching a soft gluon to the c quark line in Fig. |l|a. 
This matrix element is of relative order v 4 ~ 5%. In NRQCD, matrix elements such as 
(Q\xlD 2 tp c \B c ) and (0| xlgff ■ Bip c \B c ) are universal quantities which can, in principle, be 
extracted from other B c decays. Note that the relativistic contributions are formally smaller 
than QCD radiative corrections of order a s , which we have not included.! 



IV. VERY SOFT PHOTONS AND THE B* c POLE 

For weak radiative B c decays, there is the possibility that the outgoing photon is very 
soft, so much so that there is a large time separation between the event where the photon is 
emitted and the event where quarks annihilate to the lepton pair. In this case, the physics is 
not adequately described by (even nonrelativistic) quark fields, and consideration of hadronic 
intermediate states is necessary. Let us study this part of the spectrum more carefully. The 



treatment in this section is similar to that of Ref. [12| for the decay B — > ixlv. 

The amplitude A(B C — > Ivy) is a second-order process involving contributions both 
from the electromagnetic part of the Hamiltonian H em and from the weak part H w . For 
a sufficiently soft photon, in which case the recoil of the hadrons can be neglected, the 
amplitude is 



A(B C -> Iwy) = % J d'xe 1 ^ (lu^r\T{H el (x),H w (0)}\B c ) 



= Y,{lv\H w \M) 1 -^—{M 1 \H el \B c ) + local terms, (4.1) 

where AEm is the energy difference between the B c and the M7 intermediate state, and 
M = B c , B*, B**, . . . represents any meson of the B c family. Heavy quark spin symmetry 
implies that the mesons M are found in degenerate pairs, of which the B c and B* comprise 



'The leading radiative correction to B c — > Iv was computed in Ref. [11] 
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the lightest. Heavy quark symmetry also implies that in the limit E n — > 0, the matrix element 
(M^\H e i\B c ) vanishes except when M is a member of the ground state doublet (B C ,B*). 
In this limit, this is a "zero- recoil" transition, such as at the w — 1 point in semileptonic 
B decays |13 |. Therefore, for very small E 1 the sum in Eq. ( |4.1|) is dominated by the 



contributions from B c and B*. Since the weak B c — > Iv transition is helicity suppressed for 
mi = 0, only the vector B* state survives the sum. This situation is represented in Fig. ^a. 

We will neglect contributions from the graphs of Fig. |]b, corresponding to intermediate 
ip states. These graphs are proportional to g^/m^ and are formally suppressed compared 
to the B* pole. Moreover, since the ip states are far off shell, these contributions may be 
systematically accounted for by including a set of higher-order NRQCD operators. Note that 
the situation is completely different in B u — > lwy decays, where the analogous contribution 
is from an intermediate p meson and cannot be calculated in a model independent fashion 
in HQET, thus providing an intrinsic uncertainty. There are also local terms arising from 
the graphs like Fig. |]c. Concentrating on the "pole" part, 

A(B C - W) =Y.^\H W \BD- l -——(B* cl \H el \B c ) + .... (4.2) 

m ^rn Bc v ■ \—k) 

Clearly, the position of the pole in Eq. (|4.2j ), although consistent with a heavy quark sym- 
metry expectations, is not correctly located at = 0. We now show how the leading 
chromomagnetic corrections move the pole to the right place. 

At this stage it is convenient to adopt an HQET notation, in which the NRQCD operators 
are build from four- component spinors constrained by a set of on-shell conditions, 

e = %, %f = -* b . (4.3) 

In this formalism, the electromagnetic part of the Hamiltonian contains two parts, which 
can be identified with the spin-conserving electric and spin-flipping magnetic transitions, 



H em = ^ c a^F a ^ c + ^ b a^F a s^ 



4m c - -c**^4^*V> (4-4) 

where F a p = d a Ap—dpA a is the electromagnetic field strength tensor. Clearly, only the mag- 
netic part of H em will contribute to the B* intermediate state in Eq. ( |4.2| ). The {B* c ^\H em \B c ) 
matrix element is then 
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i{B*^\H ern \B c ) = -k a e* 



' Qc (B* c \V c a a ^ c \B c ) + ^k(B* c \%<j a ^ h \B c ) 



(4.5) 



2m c ' 2m& 

These matrix elements can be calculated using the trace formalism ( |3.7| ), yielding 

(B* c \%<t°0%\B c ) = 2m BcV ; Vu e^ , (4.6) 

where r] n is the polarization vector of the B* meson, and ^ q represents either heavy effective 
field. 



We now use heavy quark symmetry (|3.5| ) to relate the decay constant of the B* to that 
of the B c , 

(0|f 67/ ,(l - l^c\Bl) = if Bc m BcV , . (4.7) 
Thus in the heavy quark limit and with mi = we obtain for the soft photon amplitude 

A(B C - lisj) = -Y^ff B m Bc E*l e^l^elkp + local term, (4.8) 

V 2 V • K 

where we have introduced for convenience the perturbative B c magnetic moment, fi Bc = 
eQb/2mb + eQ c /2m c . Since in the heavy quark limit the B c and B* are degenerate, the pole 
is at the "wrong" position E 1 — 0. Next we show how inclusion of the l/m C) & corrections 
removes the degeneracy and shifts the pole to E 1 = —A = -(mg. -mgj, in the unphysical 
region where it belongs. 

The leading l/m c ^ corrections come from the insertion of the kinetic energy and chro- 
momagnetic dipole operators, Okin = ^ q (iD) 2 ^ q /2m and O mag = g s ^/ q a al3 G a/ 3^/ q /4:m. We 
use the customary definitions 

{B^\%(iD) 2 %\B^) = 2m Be X 1 
(B^\^ q a^G a ^ q \B^) = 2m Bc d$\ 2 , (4.9) 

where du = 3 and d* M = —1. Note that unlike in HQET, in NRQCD the matrix elements 
of these operators have different velocity powers; the first operator is there at leading order, 
while the second is suppressed by a power of v. 

Inserting Okin and O mag vertices on the B* meson line gives rise to the double pole 
contribution 

A(B c ^lvj) = (lu\H w \B^){B^\O kin + O mag \B^){B^ 1 \H d \B c ) ( 1 ) 

\2m Bc v ■ k J 

~ ^(lp\H w \B^)(Bi*^\H el \B c ) + ... , (4.10) 



2m Bc (v ■ k) 

where 5* = — (Ai + d* M \<2)j2m re d and m re d = mbm c /{mb + m c ). The same matrix elements 
contribute to l/m Cy b corrections to the physical meson masses, 

+ m c --^—{X 1 + d% ) \ 2 ). (4.11) 



^red 



>''■ ""' 2m 
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(b) 



FIG. 5. Pole structure in B c — * Ivy, at next-to-leading order. The Oj represent insertions of 
the kinetic energy and chromomagnetic operators. 



When Oki n and O mag and inserted on the external line as shown in Fig. f|b, they modify the 
heavy meson propagator by shifting the position of the pole, 
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v ■ k v ■ k + 5 v ■ k 



v ■ k 



+ ... 



(4.12) 



where 5 = — {\\ + du^) / ^^red- Taking this modification into account in (4.8) and combinig 
it with the double pole contribution ( 4.10|) , one sees that the terms proportional to Ai cancel 
but terms proportional to A 2 don't, since g?m ^ d* M . Then the propagator in (|4.2|) is replaced 
by 



v ■ (—k) v ■ (—k) 



1 + 



v(-k) 



v ■ k + A 



Hence we obtain 



i/n ; \ G F V cb fl B , 
A{B C -> Ivy) = —f Bc m B , 



V2 



v- k + A 



e^l„v a e* v k p + local, 



(4.13) 



(4.14) 



with the pole now in the correct place. Since the photon energy spectrum is finite, and in 
fact vanishes, as E 1 — > 0, the shift in the position of the pole has little effect on the total 
rate. 



V. SUMMARY 



We have performed a model-independent study of the weak radiative B c decay in the 
framework of nonrelativistic QCD. This process is an important competitor to the annihi- 
lation process B c —>■ \w, which eventually could be used to extract the B c decay constant. 
We found that the branching ratio for B c — > 7/i^ M is of the same order of magnitude as 
the corresponding purely leptonic decay, while B c ^yev e completely dominates its leptonic 
counterpart. We have generalized the NRQCD Lagrangian by introducing external sources 
for the electromagnetic and weak fields. At leading order, NRQCD gives a result similar 
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to what one might expect from a constituent quark model; at higher order, we estimate 
that the leading nonperturbative corrections to T(B C — > ^Iv) /T(B C — > Iv) are at the level of 
5%. Finally, we showed how the B* pole emerges for very soft photons. Unfortunately, the 
branching ratio for this process is small, B(B C — > Z1/7) ~ 4.4 x 10" 5 , and it will certainly be 
a challenge to observe this process even at future hadronic B Factories. 
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